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ABSTRACT 

Using a large suite of galaxies with a variety of concentrations and metallicity gradients, we study the evolution of non-dissipative 
("dry") equal mass mergers. Our purpose in generating these simulations is to explore how the metallicity gradient in dry mergers 
depends on the structure and metallicity gradients of the galaxies involved in the merger Specifically, we would like to answer: 
Could dry mergers lead to metallicity gradients as observed in elliptical galaxies in the local Universe? Do dry mergers always 
lead to a flattening of the initial (i.e., pre-merger) gradient? From this modeling, we conclude that: The ratio of the remnant and 
the initial galaxy slopes span a wide range of values, up to values greater than 1 (with values greater than one resulting only when 
companions have gradients twice the progenitor). For a merger between two ellipticals having identical initial metallicity slopes (i.e., 
equal companion and galaxy slopes), the metallicity profile of the remnant flattens, with a final gradient about 0.6 times the initial 
one. Ellipticals can maintain their original pre-merger metallicity gradient if the companion slope is sufficiently steep. The amount of 
flattening neither depends on the characteristics of the orbit of the progenitors or on their initial concentration. Given the diversity in 
outcomes of the mergers, we conclude that dry mergers do not violate any observational constraints on the systematic characteristics 
of metallicity gradients in local ellipticals. In fact, dry mergers may be important within the context of the results of our simulations 
and may explain the large scatter in gradients for massive ellipticals and the relative paucity of massive ellipticals with no or shallow 
metallicity gradients. 
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] 1. Introduction 

, Recently, we have seen a revision in our understanding of 

■ the formation of elliptical galaxies which was motivated by 
' the need to reconcile the apparently "anti-hierarchical" behav- 
. ior of Active Galact ic Nuclei (e.g., ,Cattan eo & Bernardi 2003;: 
' iHasinger et al.ll2005l), the evolution of their luminosity function 
I with redshift (e.g., iBundv et al.l l2006l) . as well as the various 
. lines of arguments based on the analysis of their stellar popu- 
' lations (iNelan et al.ll2005h . In particular, recent results suggest 
[ that we need a substantial modification in how we treat the rela- 

■ tionship between the baryons and the dark matter in early type 

■ galaxies. Pointedly, more massive ellipticals are older and have 
forme d faster than their lower mass counterparts ( iThomas et al.l 
l2005h . In addition, independent evidence that supports such 
downsizing is inferred from the systematic characteristics of 
abundance ratios in elliptical galaxies, namely, the increase of 
the mean [M g/Fe] in the st ellar populations of el lipticals with 
galaxy mass ( iWorthev et al.l 1992; Matteu c^l 19941) . 

Hierarchical modeling, in its most recent incarnations, only 
partly accounts for the downsizing observed in elliptical galax- 
ies. In practice in such models, significant mass assembly still 
occurs at late times, but most of the stars have been formed at 
high redshift in small subunits that merge to grow early type 
galaxies. The currently preferred mechanism for the assembly 
of massive spheroids is a sequence of "dry" merger (e.g., 
iKhochfar & Burkertll200llDe Lucia etai]|2006l: ICattaneo et al.l 



' In this paper dry merger means a pure dissipationless merger of an 
elliptical stellar systems, i.e., without any gas or significant subsequent 
star formation. 



120061 12008h . There is evidence that dry mergers might solve 
some of the outstanding issues related to the growth and evo- 
lution of elliptical galaxies. For instance, the brightest ellipti- 
cals have central phase-space densities comparable to those of 
disk galaxies, suggesting that, if mergers are responsible of their 
build-up, they likely req uire only a small amount, if any, dis- 
sipation (ICarlberdll986l) . Moreover, dry mergers might explain 
the strong size evolution of massive red galaxies with redshift 
(see, e.g., Ivan Dokkum et al. 2008; Khochfar & Silk 2006) and 
the formation of s lowly rotating ellipticals with boxy isophotes 
(iNaab et al.ll2006l) . 

One important aspect of determining whether or not dry 
mergers are a viable mechanism for explaining the evolution 
of elliptical galaxies is how various aspects concerning the 
metal content and d istribu tion changes after a merger event. 
iPipino & Matteuccil (l2008l) show that a series of multiple dry 
mergers (with no associated star-formation) involving building- 
blocks that satisfy the [Mg/Fe]-mass relation c annot fit the mass- 
metallicity relation and vice-versa (see also ICiotti et akl Il99ll; 
ICimatti et al.ll2006l ). A major dry merger (mass ratios approxi- 
mately about one), on the contrary, does not violate these ob- 
servational constraints if such a merger occurs between galaxies 
that already obey both the mass(cr)-[Mg/Fe] and the mass(cr)- 
metallicity relations. However, because this mechanism would 
only operate over a limited range of initial masses, meaning, you 
cannot have a long series of major mergers to make the most 
massive ellipticals that initially obeyed the mass-metallicity re- 
lation, thus this process alone cannot be the underlying physical 
cause of these trends. 
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But global trends and relationships amongst ellipticals is 
not the only constraint available in testing the validity of the 
dry merger scenario - the structure within individual galax- 
ies can provide further significant constraints. Radial negative 
metallicity gradients a re a common feature in the stellar pop - 
ulations of spheroids dCaroUo et al.lll993l: iDavies et anil993l) . 
Observations show that the majority of ellipticals have a typi- 
cal decrease in metal Ucity of 0.3 dex per decade in radius but 
with a large scatter dAnnibali et al.l l2007l) . Moreover, a posi- 
tive correlation between the metallicity gr adient with th e galac - 
tic mass has been originally claimed by ICaroUo et alj (Il993l) . 
but only for masses lower than ~ IO^Mq. More recently, addi- 
tional studies have also found a positive co rrelation of the gra- 
dient with rnass for a wide range of masses dOgando et al.ll2005l: 
Forbes et al.l l2Q05^.'Sanchez-Blazque z et alJl2007l:ISpolaor et al. I 



20091 but seeLKobayashi & Arimoto 1999HAnnibaU et al.ll2007l) " 



Steep metallicity gradients are expected f r om classical dis - 
sipative collapse models (e.g., iLarsonl 119741: ICarlberd Il985l) . 
More modern versions of dissipative collapse starting from 
(semi-) c osmological boundary/initial conditions dKawatal 19991: 
|Ch iosi & Carraro 2002) can also explain the relation between 
metallicity gradient slope and stellar mass. The predicted metal- 
licity gradient can be as high as -0.5 dex decade ' in radius, 
hence galaxies can exhibit a slope steeper than the average. 
Owing to such a steep predicted slope it has been hypothesised 
that galaxy mergers and monolithic collapse c ould be at work to- 
gether in order to explain the range in slopes (iKobav ashi 2003). 
On the other hand, more realistic models with a detailed treat- 
ment of the chemical evoluti on predict a typica l metallicity gra- 
dient of -0.25 dex decade ' (iPipino et al.ll2008b and also a trend 
between the gradient and mass. These ar e just the trends n ec- 
essary to explain the latest observations dPipino et al.ll2009l in 
preparation). This finding implies that the assumption that ellip- 
ticals are formed through a mixture se veral channels (inc luding 
mergers and monolithic collapse, e.g. iKobava shT 2004) is not 
necessary, but cannot be excluded with current available data. 

N-body numerical simulations can be an important tool for 
understanding the evolution of metallicity gradients in mergers 
and their remnants. Such simulations are able to account for dif- 
ferent environments, mass ratios, and morphologies of the pro- 
genito r galaxies. The cosmological simulations by iKobavashil 
d2004l) have shown that a variety of gradients in early-type galax- 
ies can be generated and are mainly due to a difference in merg- 
ing histories. In such a scenario, gradients are destroyed during 
mergers, by an amount depending on the mass ratio of the pro- 
genitors. Gradients can be regenerated if strong star formation 
occurs in the central regions - thus requiring sufficient gas dis- 
sipation to fuel such star-formation - and they slowly evolve 
through subsequent gas accretion. These simulations indicate 
that, in order to reproduce the large scatter observed in metal- 
licity gradients, both gas infall (monolithic collapse) and major 
mergers must occur While cosmological simulations can follow 
a number of different merging histories (from the coalescence of 
small gas-rich subunits to major dry or wet mergers), it is also 
important to understand the role that each of these processes play 
in reshaping any initial metallicity profile. 

While a limited number of studies have investigated 
the evolution of metallicity gradients in gas-r ich mergers 
dMihos & Hernquisllll994t iBekki & Shiovalll999l) . taking into 
account both dissipative processes and the impact of star forma- 
tion, less attention has been devoted to the evolution of metallic- 
ity profiles in dry mergers. Early attempts at these types of sim- 
ulations suggested that metallicity gradients flatten, due t o both 
mixing and core-halo differentiation during the merger dWhitd 



Il980h . However, N-body experiments have shown that the ini- 
tial state of the galaxy is not completely washed out during the 
coalescence of the systems (van Albada 1982), suggesting that 
radial abundance s may be only moderately reduced by d issipa- 
tionless mergers dBarneslll996l:lMihos & Hernquistlll994 . 

Given the importance of this problem for understanding the 
evolution of early type galaxies, we wish to investigate impact 
of dry mergers on galaxy evolution further through the use of 
N-body simulations of equal-mass elliptical galaxies. From this 
investigation, we want to understand: 1) if dry mergers always 
cause the flattening of an initial galaxy gradient; 2) under what 
conditions does this flattening occur; 3) what is the magnitude 
of the flattening that is typically produced for a range of initial 
conditions; and 4) how does this evolution depend on orbital and 
morphological properties of the interacting systems. To achieve 
this aim, we simulated a number of equal-mass dry mergers, with 
a variety of orbital parameters, morphologies, and initial metal- 
licity gradients for each component of the merger. We have cho- 
sen to begin investigating equal mass mergers, because they are 
likely to result in the greatest amount of mixing in phase-space 
and so, ultimately, they should cause the strongest variations in 
the initial (i.e. pre-merger) metallicity profiles. 

The main characteristic of the models (initial conditions and 
numerical code used) are briefly described in Sections |2] and [3] 
results are presented and discussed in Section|4]and conclusions 
are drawn in Section |5] 



2. The model 

The mergers studied in this paper consist of the coalescence of 
two equal-mass elliptical galaxies, without any dissipative com- 
ponent. The fiducial case is described in Sect. 14. 11 It corresponds 
to mergers of two giant elliptical EO (hereafter gEO) galaxies, 
whose internal parameters are given in Table[T] moving on 6 dif- 
ferent orbits (2 parabolic, 1 hyperbolic, 1 radial and 2 ellipticals; 
see Table |2]i. 

To understand to what extent the evolution of the metallic- 
ity gradients depends on the galaxy parameters adopted, we also 
modeled interacting galaxies having a central density 5 times 
higher and 2.5 times lower than that of the gEO galaxies (we refer 
to these systems as gEOm and gEOl, respectively; see Table[T]for 
their morphological parameters). To allow the models to reach an 
equilibrium and to demonstrate their initial stability, all the mod- 
eled galaxies were allowed to evolve in isolation for 1 Gyr be- 
fore any encounter This allows us to distinguish internal relax- 
ation phenomena from those due to the merging process itself. 
These systems represent quite idealized initial conditions: they 
are completely spherical, showing no rotation in their baryonic 
component. Observations show, however, a vas t range of ellip- 
ticities and levels of rotation (see, for example lEmsellem et al.l 
|2007|) . The systematics of the amount of rotation and elliptic- 
ity with mass seem to suggest that more massi ve spheroids are 
not likely to be the resul t of dissipative mergers dCarlberdl 19861: 
Ivan Dokkum et al]|2008l) . Therefore, to test "more realistic" ini- 
tial conditions, we modeled also mergers between two equal- 
mass galaxies (hereafter called iEO), having masses equal to 
0.5 times that of gEOs. These galaxies were then subsequently 
merged to form a system (hereafter called 2iE0) whose total 
mass is equal to that of gEOs, but having initially a diff'erent mor- 
phology (triaxial instead of spherically symmetric) and a certain 
amount of rotation (see Table [1] for the initial parameters of the 
iEO galaxies and Fig.[T]for some properties of the initial galaxy 
models). In this way, we were able to, at least partially, account 
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Fig. 1. Initial properties of some galaxy models: (a) projected-density map of the gEO galaxy, after evolving in isolation for 1 Gyr 
Contours are equally spaced on a logarithmic scale; (b) projected-density map of the 2iE0 galaxy which is the result of a collision 
between two iEOs elliptical models (see text for details). We used the same scale interval and contours level as in the previous 
panel; (c) volume-density profiles of some galaxy models, as indicated in the legend of the figure; (d) line-of-sight velocity and 
velocity dispersion profiles for the gEO galaxy model (black solid and black dashed lines respectively) and for the 2iE0 galaxy 
merger remnant (gray solid and gray dashed lines respectively). 



for more complicated merger histories and hopefully more real- 
istic merger components in our simulations. 

Each galaxy in the merging pair has been modeled with 
- 120000 particles, distributed among stars and dark mat- 
ter. We also re-simulated some cases, using a number of par- 
ticles two and four times higher, attempting to insure that the 
results presented here do not depend on the numerical resolution 
adopted in our models. 

Finally, for each of the galaxies involved in the interaction 
and merger, we have adopted 15 different initial metallicity pro- 
files. These profile were constructed by assigning a metallicity, 
Z, to each star particle in the system that depended exponentially 
on the distance of the star from the galaxy center (the formal 
functions are given in TableO. This is also why the galaxies are 
allowed to evolve in isolation for 1 Gyr This time allows us to 
ensure that the metallicity gradient of the isolated galaxy is sta- 
ble and remains as it was assigned. A large variety of initial cen- 
tral metallicities and gradient slopes are thus studied in this way, 
ranging from flat profiles (as it is the case of gradient id=grad02 
and grad05) to metallicity gradients around -0.3 dex per decade 
in radius (gradOl, gradll, gradl2) to quite steep slopes (gradient 
id=gradl4, gradlS). In this and in the next sections, the metal- 
licity gradient is defined as. 



A/og(Z) 
Alog(r) 
log{Z{Q.lr,o))- 



(1) 



■ logiZir^o)) 



logiZ{rso))-log(z(0.lrso)), 



where rso is defined as the radius which contains half of the 
baryonic mass ("baryonic half-mass radius"). 
In total, we analyzed 28 diff'erent simulations of major dry merg- 
ers, and, for each interaction, we studied 120 different combi- 
nations of the initial metallicity profiles, corresponding to 3360 
different possibilities. 

3. Numerical method 

All the sim ulations have been run usi ng the Tree-SPH code 
described in 'Semelin & Combes" (2002!). A further description 
is give n also in Di Matteo et al. ( 2007) and in Di Matteo et a L 
(|2008|) . The code uses a hierarchical tree method fearnes & Hut 
Il986l) to evaluate gravitational forces and a smoothed particle 



Table 1. Galaxy parameters. Both the stellar and dark mater pro- 
files are represented by Plummer models, having characteristic 
masses, respectively, given by Mb and Mh and core radii given 
by rs and rn 





gEO 


gEOl 


gEOm 


iEO 


Mb 


[2.3 X 10'' Mo] 70 


70 


70 


35 


Mh 


[2.3 X 10' Mo] 30 


30 


30 


15 


re [kpc] 4.0 


6.0 


2.0 


2.8 


rH [kpc] 7.0 


7.0 


7.0 


5.0 


N„ar 80000 


80000 


80000 


40000 


Ndm 40000 


40000 


40000 


20000 


Table 2. Orbital parameters 


orbit id 


ini ^'iiii 


L" 








[kpc] [100 kms-^] 


[lO-kms-' 


kpc] 


[10^ km^s-^] 


01 


100. 2.0 


57.0 




0. 


02 


100. 3.0 


59.0. 




2.5 


05 


100. 2.0 


80.0 




0. 


13 


100. 0.0 


0.0 




-2. 


14 


70.0 1.5 


41.3 




-0.66 


15 


70.0 1.6 


41.3 




-1.57 


" It is the absolute value of the anj 


^ular momentum of the unit mass 


i.e., L 


=1 r X V |. 









It is the total energy of the relative motion, i.e., 

E = v'/2 - G{m, + m2)/r, with m, = m2 = 2.3 X 10"Mo 



hydrodynamics method to s imulate the evolution of gas iLucyl 
(ll977h : lGingold & Monaghan (1982). Since the work presented 
here only investigates dry-mergers, only the part of the code 
evaluating the gravitational forces acting on the systems has 
been used. 

Gravitational forces are calculated using a tolerance param- 
eter 6 - 0.7 and includes orders up to the quadrupole term in 
a multiple expansion. A Plummer potential is used to soften 
the gravity at small scales, with constant softening lengths of 
e = 280 pc for all particles. The equations of motion are in- 
tegrated using a leapfrog algorithm with a fixed time step of 
0.5 Myr. 
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Fig.2. Evolution of the metallicity profiles in some major (1:1) dry mergers of two gEO galaxies. Metallicity is shown versus 
the distance from the galaxy center for the two galaxies before the interaction (black and thick gray solid lines) and for the final 
remnant for different orbits (see figure legend). The values on the abscissa are shown in units of the half-mass radius r^o of the 
corresponding galaxy. Different panels in the figure correspond to different initial metallicity profiles for the two interacting galaxies; 
(top-left panel) id=grad01grad01; (top-middle panel) id=grad01grad02; (top-right panel) id=grad01grad03; (bottom-left panel) 
id=grad01grad04; (bottom-middle panel) id=grad01gradl5; (bottom-right panel) id=gradl5gradl5. 



Table 3. Initial gradients and metalUcity profiles for the elUptical 
galaxies 



Gradient id 


A 


Metallicity profile 


gradOl 


-0.26 


Z(r) 




3zolO-"-'"^ 


grad02 


0.00 


z{r) 




Zq 


grad03 


-0.13 


z{r) 




3zolO-°"35. 


grad04 


-0.13 


z{r) 




1.5zolO-""35 


grad05 


0.00 


Z{r) 




1.5zo 


grad06 


-0.22 


z{r) 




3z0lO-o"«'- 


grad07 


-0.18 


z{r) 




3zolO-''°5'- 


gradOS 


-0.15 


z{r) 




3zolO-°"*'" 


grad09 


-0.07 


zir) 




3zolO-°''2^ 


gradlO 


-0.04 


z{r) 




3zolO-»-'"^ 


gradl 1 


-0.30 


z{r) 




3zolO-»-''**^ 


gradl2 


-0.33 


zir) 




3zolO-»-'»^ 


gradl 3 


-0.37 


zir) 




3zolO-o'o^ 


gradl 4 


-0.41 


zir) 




3zolO-o"' 


gradl 5 


-0.44 


zir) 




3zolO-o'2^ 




Initial slope 

Fig. 3. Absolute values of the initial and final metallicity slopes 
for gEO-gEO mergers. Different orbits are indicated by different 
symbols, as explained in the legend. The initial slope has been 
evaluated by superposing the initial metallicity profiles of the 
two progenitor galaxies. 



4. Results and discussion 

4.1. Evolution of metallicity profiles in major dry mergers 

The evolution of some of the metallicity profiles during major 
non-dissipative mergers are shown in Fig|2] In this plot the 
initial morphologies of the two interacting galaxies are always 
the same, while the orbits and initial metallicity gradients have 
been varied. In some cases, as for the top-left and bottom-right 
panels, the two interacting galaxies have initially the same pro- 
file. This situation of having identical initial equal metallicity 
profiles is the case studied by White (1980): in this cases, the 
profile flattens, due to the redistribution of stars in the systems 
induced by the merging process. Indeed, one can clearly see 



from the final profile that the flattening is due to both the 
decrease of the central metallicity and to the contemporaneous 
increase of the metallicity in the outer regions. Note also that 
the evolution of the profile does not seem to be very sensitive to 
the orbital initial conditions: orbits with different initial energies 
and angular momenta produce very similar final profiles, at 
least in the region inside the remnant half-mass radius. This 
happens because mixing and phase-space density evolution are 
independent on the orbital parameters, as discussed in Appendix 

m 

Fig.|2]also shows some examples of the two galaxies having 
diff'erent initial metallicity gradients and distributions ("mixed"). 
As discussed in the Introduction, observations (see, for example 
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['•merging] M/M^g^ M/M^g^ 



Fig. 4. Core-halo differentiation and mixing in some of major gEO-gEO mergers modeled for this study, (left) Shows the evolution 
with time of the Lagrangian radii containing 10%, 25%, 50%, 75%, and 90% of the total stellar mass. For simplicity, the radii are 
shown only for stars belonging to one of the two galaxies (the behavior is the same for the companion). The time axis is in units 
relative to the merging time. Mergers with a variety of orbital characteristics are shown, as indicated into the figure legend. {Middle 
and right) Metallicity profiles for some gEO-gEO mergers, as a function of the fractional mass contained. The panels correspond 
to different metallicity profiles: id=grad01grad01 (middle panel); id=gradl5gradl5 (right panel). Metallicity is shown for the two 
galaxies before the interaction (black and thick gray solid lines) and for the final remnant. 



lOgando et all 120051: [Spolaor et al. Il2009[) show that early-type 
galaxies in the local Universe have a wide dispersion in their 
metallicity gradients, particularly at the high end of the distribu- 
tion of masses. This means that dry mergers between galaxies 
with different metallicity gradients should be common locally. 
Of course, this is not true if there is some sort of regularity in 
or correlation between the metallicity gradients as a function of 
loca l galaxy density, for which yet there is no clear e vidence 
(see lslnchez-Blazquez et"aLll2006l: IClemens et al.ll2009b . But if 
"mixed" mergers are common (as seems likely), then it is in- 
teresting to understand what the resulting gradient is and how 
it depends on the initial profiles and orbital parameters of the 
progenitors. By examining the plots (Fig. |2]i and analyzing the 
simulations, the resulting gradients in mixed mergers can lie in 
between those of the progenitors: the merging of galaxies with a 
flat and steep gradient (A = -0.26) produce a remnant having a 
slope of about -0.09 which is lower than the steepest progenitor 
gradient but certainly not flat as in the companion galaxy. This 
implies that for a given progenitor galaxy with an initial metal- 
licity gradient, major dry mergers do not necessarily lead to a 
flattening of the initial galaxy slope, but the outcome depends 
greatly on the slope of the companion galaxy. An early-type el- 
liptical galaxy merging with an other having the same initial 
metallicity profile will have a final slope lower than the initial 
one - the amount of flattening is evaluated in Sect j4.31 - but if it 
merges with a companion having a sufficiently steep profile, the 
final slope can be the same or steeper than the initial one. 
Comparing the superposition of the initial metallicity profiles of 
tiie two progenitor galaxies (i.e. Z\(r) + Z2(r)) with the remnant 
profile, we find that, on average, the final metallicity gradient 
is a factor 0.6-0.7 flatter than the initial one, with only a small 
dispersion around the average (see Fig|3]l. 



4.2. Core-halo differentiation and phase-space mixing 

In Sect l4.1l we showed that the metallicity gradient of the final 
merger remnant depends on the initial gradients of the progenitor 
galaxies. Here we present more details about the mechanisms 
causing the evolution of these profiles. Importantly, we note that 
the metallicity profiles given in Fig |2] are shown in units of the 
stellar half mass radii r^o of the progenitors and remnant. Thus, 
a variation in the half-mass radii of the remnant with respect 
to that of the progenitor galaxies, without any variation in the 



Z(r) profile, should be sufficient, in principle, to change the final 
metallicity profile. 



Indeed, as shown for gEO-gEO mergers in the left panel of 
FiglH the half-mass radius changes with time, at each pericenter 
passage and mostly during the final phases of coalescence of the 
two systems (t/tmergmg -1), with the final value being about 50% 
greater than the pre-interaction value. Note also that outer parts 
of the systems (those outside r^o) show the strongest evolution. 
Indeed, the tidal shocks increase the radius rgo, containing 90% 
of the stellar mass by more than a factor of 3. On the contrary, the 
inner radii show only modest change after the merger Moreover, 
the evolution with radius also has a behavior that depends on the 
orbital geometry of the progenitors. The dynamical evolution of 
the outer radii is affected by the orbital geometry while that of 
the inner regions is largely independent of the orbital parameters. 



As first shown bv I White! (Il980h . it is possible to separate the 
effects due to structural changes from those arising from mixing 
by showing the metallicity profiles as a function of the enclosed 
mass rather than as a function of some physical coordinates. If 
the evolution of the metallicity gradient was only due to core- 
halo differentiation, then the initial and final profiles in enclosed 
mass would be identical. However, a clear flattening occurs dur- 
ing the merger (see Fig. |4] middle and right panels), indicating 
that mixing must contribute to the evolution of the metallicity 
gradient. In our si mulations, thi s flattening is more pronounced 
than that found bv lWhitel(ll980l) . This difference could be due to 
his choice of concentrating on merger remnants which showed 
density profiles and velocity structures as close as pos sible to 
that of the progenitors (see discussion in IWhitd [19801) . In ad- 
dition, the process of mixing does not seem to depend on the 
initial orbital geometry but it is generic to the merging process 
itself and leads to rather similar metallicity gradients. This im- 
plies that mixing is mostly driven by the progenitor struct ure that 
(at least partially) survives the merging process (see also lWhit^ 
I1980H b arnes l ll996HValluri" et al.ll2007l) . As a consequence, our 
simulations indicate that it is not possible to disentangle the or- 
bits of the parent galaxies by using the metallicity distribution of 
the remnant if the analysis is restricted to the central region of 
the galaxy (i.e., inside the half-mass radius). 
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Fig. 5. Remnant (absolute) slopes as a function of the initial 
galaxy (absolute) slopes, for different orbits and different ini- 
tial galaxy models. Symbols '01', '02', '05', ' 13', ' 14' and '15' 
in the legend correspond to interactions among two gEO galax- 
ies moving on different orbits (see Table[2]for the corresponding 
orbital parameters); 'Oil' and '021' correspond to interactions 
among two gEOl galaxies with orbit id=01 and 02, respectively; 
'01m' and '02m' correspond to interactions between two gEOm 
galaxies with orbit id=01 and 02, respectively. Elliptical galaxies 
resulting from the merger of two 2iE0 galaxies are shown (with 
the symbols '01,2iE02iE0' and '01,2iE02iE0' referring to or- 
bit id=01 and 02, respectively). For comparison, also the slope 
resulting from the merger of two iEO galaxies - having a total 
mass half of that of gEO - are shown. Error bars correspond to 
the standard deviation of the mean. 
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Fig. 6. Left panel: Evolution of the metallicity profiles for three 
different major dry mergers, involving respectively two gEO 
galaxies (black solid line), two gEOl galaxies (dashed red line) 
and two gEOm galaxies (dotted blue line). The initial profile is 
shown with corresponding thin curves. All the profiles are shown 
as a function of the distance r form the galaxy center. Right 
panel: Same as the previous panel, but this time all the profiles 
are shown in units of the half mass radius r^o of the correspond- 
ing galaxy. 



4.3. Metallicity gradients in major dry mergers: flattening 
steepening? 



or 



After describing how metallicity profiles evolve during a merger 
and the physical processes which determine this evolution, we 
now want to quantify the change in the radial profiles themselves 
as a function of the initial metallicity distribution. 

In Table H) we tabulate the initial and final metallicity 
gradients for simulations of two gEO galaxies with a specific 
orbital geometry (id=01) but for a range of metallicity gradients. 
To generalize this, Fig|5]shows the absolute values of the initial 
and final gradients from the simulations for all the galaxy 
morphologies and orbits. For a given initial metallicity gradient 
there are a range of final gradients, depending mainly on the 
profile of the companion galaxy. The initial galaxy models span 
a wide range in metallicity slopes ,0.0- -0.7, while the remnant 
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Fig. 7. Ratio of the companion (initial) slope to the galaxy (ini- 
tial) slope versus the ratio of the remnant slope to the galaxy 
(initial) one. Symbols are described in caption of Fig|5] Error 
bars correspond to the standard deviation of the mean. Note 
that the galaxies having a remnant slope of zero are not plot- 
ted in this figure. This corresponds to the situtation where the 
initial galaxy slope and companion slope are both zero and thus 
l/(galaxy slope) is infinite. 



slope only spans about 0.0- -0.4. Also, a close inspection of 
Table |4] tells us that the final slope is always flatter than the 
maximum between the parent and his companion slopes. In most 
of the cases, the final slope is lower than the mimimum between 
the parent and his companion slopes. Note that, as suggested 
previously, galaxies having the same morphologies have very 
similar final slopes, independent of their orbital geometry. This 
is the case regardless of the initial morphologies. In addition, 
the most concentrated galaxies (gEOm) have the smallest range 
of gradients and generally smaller gradients for the same initial 
conditions while galaxies with low concentrations (e.g., gEOl) 
generally have larger slopes and exhibit a wider range of slopes. 
This is due to the fact that we have assigned to each galaxy the 
metallicity profiles following a simple Z=Z(r) law where r is the 
physical distance from the galaxy center of mass. That is, the 
initial metallicity profiles have not been normalized in anyway 
to the galaxy half-mass radii, while the metallicity gradients, 
according to Eq[T] are evaluated using them. 
A comparison of the evolution of the metallicity profiles for 
galaxies with different initial concentration is shown in Fig|6] 
In the left panel, the initial metallicity profile is shown as a 
function of r, the distance from the center, for the gEOl, gEO and 
gEOm galaxy. This initial profile is obviously the same for the 
three galaxies, accordingly to our choice. The final profile of the 
remnant galaxy, in turn, depends on the initial galaxy density, in 
the sense that the more concentrated the galaxy is, the higher is 
the value of Z(r), for any given r. This is simply a consequence 
of the fact that the more concentrated the galaxy is, the more 
metal rich it is too, according to the way the initial metallicity 
profile has been assigned. When normalizing to the half-mass 
radius (right panel in Fig|6]l, in turn, the initial profiles are 
obviously different, as well as the final ones, but interestingly 
the amount of flattening of the final gradients with respect to the 
initial ones are very similar If the two galaxies have initially 
identical metallicity profiles, as in the case shown in Fig|6] the 
final metallicity gradient will have a slope equal to 0.6-0.7 times 
the initial one, with no clear correlation with the initial central 
density of the galaxy. 
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Fig |5] shows that a range of initial gradients gives rise to a 
variety of remnant slopes, this depending, of course, on the slope 
of the companion galaxy, as well as on the initial concentration. 
What is the relationship between the remnant and initial galaxy 
slopes? This is shown in Fig]?] The relation between the initial 
and remnant gradients has generally a small scatter about a mean 
relation, with higher dispersion for the high ratios. This figure 
basically summarize the results of this work which implies: 

- The ratio of the remnant and the initial galaxy slopes span a 
wide range of values, up to values greater than 1 . 

- For a merger between two ellipticals having identical initial 
metallicity slope (i.e., companion slope=galaxy slope), the 
metallicity profile of the remnant flattens, with a final gradi- 
ent about 0.6 times the initial one. 

- Remnants can have metallicity gradients greater than that of 
the progenitor elliptical galaxy. Our simulations show that 
this happens every time the companion has a slope two times 
greater than the parent. 

- Ellipticals can maintain their original pre-merger metallicity 
gradient if the companion slope is sufficiently steep. 

- The final remnant gradient does not depend on the orbits of 
the progenitor ellipticals. 

- While the final gradient does depend on the initial concentra- 
tion of the progenitors, the amount of flattening (i.e. remnant 
slope/progenitor slope) does not. 

4.4. The Role of Dissipationless Mergers in the Evolution of 
Elliptical Galaxies 

Do these results in any way suggest that equal mass dry mergers 
might be inconsistent with the metallicity gradients observed in 
local ellipticals? The short answer is no. There are no system- 
atic trends in magnitude of metallicity gradients in ellipticals 
save one, metallicity gradients generally become steeper with 
increasing mass, bu t with a large scatter (iQgando etal J 120051: 
ISpolaor et al. |2009[) . Our models show that, as long as the metal- 
licity gradient of the companion is not more than twice that of 
the parent, the remnant will have a shallower gradient. Thus one 
can imagine that the build up of the most massive ellipticals is 
through a series of dry mergers whereby the gradient gradually 
softens leading to a trend, albeit with large scatter, between mass 
and metallicity gradient. But if there are several mergers in se- 
quence with long times between then, the metallicity gradients 
would increase and decrease rather randomly, thus leading to a 
situation where there are few systematic trends between the char- 
acteristics of the metallicity gradient and structure properties of 
ellipticals. This is precisely shown in Fig (8] where the evolution 
of an initial metallicity gradient, due to recursive dry mergers, is 
followed. If dry mergers can occur among ellipticals having dif- 
ferent gradients, randomly distributed in the shaded areeQ shown 
in the plot, a sequence of such events does not lead to a clear 
trend (for example, the more massive is the galaxy, the shal- 
lower its profile). In turn, a sort of random walk, with increasing 
and decreasing slopes, seems to be the general outcome (upper 
panel in Fig|8]l, unless at each merger epoch, a given galaxy is 
allowed to merge only with companions having the identical gra- 
dient (lower panel in Fig|8j. Because of the stochastic nature of 
the process of dry mergers, it is difficult to gauge the frequency, 
and therefore the importance of this process in galaxies. While, 
most likely, the merger of two similar ellipticals would lead to 



^ This area encloses the negative metallicity slopes found by Ogando 
et al. (2005), and as shown in their Fig. 2. 



a lessening of the metallicity gradient, the final product depends 
on the concentration of the individual galaxies and their initial 
gradients. As the gradient of the parent galaxy because flatter 
and flatter, it then becomes likely that another dry merger will 
actually increase its metallicity gradient. If some of the merg- 
ers are not completely dissipationless, then this will also tend 
to increase the gradient. However, the relative mass of gas that 
forms stars during any gas accretion or merger must be small. If 
the fraction of mass forming stars were large or did so over ex- 
tended periods of time, then this process would violate the well- 
known enhancement in [a jFe] in ellipticals, especially the most 
massive ones. Therefore, any enhancement in the gradient due 
to star-formation must be small and the influence of the dissi- 
pationless component must be dominant in determining the final 
metallicity gradients in massive elliptical galaxies. 

One of the obvious outcomes of all of this is that if ellipticals 
at all masses were formed with the same gradient then we would 
expect the scatter in the metallicity gradient as a function of mass 
to increase with both time and mass. This is simply because our 
simulations show that metallicity gradients that are initially steep 
will most likely soften. This will continue up to the point where 
the gradient is shallow enough to make it more likely that the 
merger will be with a galaxy that has a slope steep enough (by a 
factor of 2) to actually increase the gradient in the remnant. This 
variation in outcomes would thereby have the eff'ect of increasing 
the scatter as galaxie s grow more rnassive . This is similar to what 
has be en observed (lOgando et alj|20()5l but see lAnnibali et al.l 
(l2007l) ). It may also explain the relative lack of galaxies with 
no gradients among the m assive ellipticals (e.g.. lOgando et all 
120051: lAnnibali et"al1l2007b . If the gradients became flat, then 
they are likely to merge with galaxies that are steep enough to 
increase their slopes. 

Note that, in this analysis, the metallicity gradients have been 
chosen to reflect those of elliptical galaxies in the local Universe. 
The progenitors of the most massive ellipticals observed today 
probably were low and intermediate mass ellipticals that merged 
at early epochs, and possibly had structural properties different 
from those of local galaxies. Unfortunately, the properties of dis- 
tant early type galaxies are not well determined. This leaves no 
recourse but to use the observed properties of local galaxies. 
Thus we caution that our analysis may not capture in detail the 
cosmological evolution of elliptical galaxies. 

5. Conclusions 

By means of a set of N-body simulations involving equal mass 
elliptical galaxies with no gas, we attempted to understand if 
dissipationless mergers always cause the flattening of the initial 
galaxy gradient, under what conditions does flattening happen, 
are there situations where a gradient could become steeper, what 
is the amount of flattening or steeping typically produced, and 
how does this evolution depend on orbital and morphological 
properties of the galaxies in these interacting systems. Our pur- 
pose is to explore how the metallicity gradients in dry mergers 
depends of the structure and metallicity gradients of the galaxies 
involved in the merger Specifically, we would like to answer: 
Could dry mergers lead to metallicity gradients of the ellipticals 
observed in the local Universe? 

From this modeling, we conclude that: 

1 . The ratio of the remnant and the initial galaxy slopes span a 
wide range of values, up to values greater than 1 (with values 
greater than one resulting only when companions have gradients 
twice the progenitor). 

2. For a merger between two ellipticals having identical initial 



8 



Di Matteo et al.: Gradients and major dry-mergers 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




I I I I I I I I I I I I I I I I I I I I I I I \ I 

9 9.5 10 10.5 11 

10g(M/Mg) 

Fig. 8. Evolution of metallicity gradients through a sequence 
of major dry mergers. The initial elliptical galaxy has 
log(M^l0)=9.3 and A=-0.05. (Upper panel) Different colors 
correspond to different merging histories, with companion gra- 
dients randomly distributed in the shaded area (cf. Fig. 2, right 
panel, in Ogando et al. 2005). (Lower panel) Different colors 
correspond to different merging histories, with companion gra- 
dients always equal to zero (black triangles), to the initial galaxy 
gradient (light gray triangles), and to the maximum slope al- 
lowed for a given mass, i.e., lower edge in the shaded area (dark 
gray triangles). 



metallicity slope (i.e., companion and galaxy slopes being the 
same), the metallicity profile of the remnant flattens, with a final 
gradient about 0.6 times the initial one. 

3. Ellipticals can maintain their original pre-merger metallicity 
gradient if the companion slope is sufficiently steep. 

4. The final remnant gradient does not depend on the orbits of 
the progenitor ellipticals. 

5. While the final gradient does depend on the initial concentra- 
tion of the progenitors, the amount of flattening (i.e. remnant 
slope/progenitor slope) does not. 

Given the diversity in outcomes of the mergers, we conclude 
that dry mergers do not violate any observational constraints on 
the systematic characteristics of metallicity gradients in local el- 
lipticals. In fact, dry mergers might also explain many of the 
characteristics of metallicity gradients. There are no systematic 
trends in magnitude of metallicity gradients in ellipticals save 
one, metallicity gradients generally become steeper with increas - 
ing mass, but with an increasing dispersion dOgando et al .I2005h . 
This observation is logically explained within the context of our 
simulations: dry mergers lead to an overall decrease in the gra- 
dient, if the companion slope is not steep enough (companion 
slope/galaxy slope less than 2). We note however, that as the 
metallicity gradient gets flatter, it becomes easier for an ellipti- 
cal to merge with another elliptical with sufficiently steep slope 
to actually increase the final metallicity gradient of the remnant. 
The reader should note that this mechanism requires the exis- 
tence at any given time of progenitors with enough steep slopes, 
perhaps formed through a pure monolithic channel. Thus, if dry 
merging is important, we would expect the variance in the metal- 
licity gradient to be the largest amongst the most massive ellipti- 
cal consistent which is what is observed. Moreover, since steep- 
ening can occur, especially for ellipticals with no or shallow gra- 
dients, such a model also predicts that ellipticals will only rarely 
show very shallow gradients. This is something that has perhaps 



been observed dOgando et alj|2005l: lAnnibali et al.ll2007h and is 
therefore an argument in favor of dry merging, not against it. 
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Fig.A.l. Evolution of the metalUcity profiles in some major 
mergers involving two 2iE0 galaxies. Metallicity is shown ver- 
sus the distance from the galaxy center for the two galaxies be- 
fore the interaction (black and thick gray solid lines) and for the 
final remnant for different orbits (see legends). The values on 
the abscissa are shown in units of the half-mass radius r^o of the 
corresponding galaxy. 




Fig.A.2. Evolution of the spherically averaged Q{r) - 
p{r)/<T^{r) during major dry mergers. Top panel: The profile 
Q{r) is shown as a function of the distance r (in kpc) from the 
galaxy center for the progenitor gEO galaxy (solid thick black 
line) and for the remnant galaxies, for different orbits (see leg- 
end). Bottom panel: Same as the previous panel, but this time the 
merger involves two 2iE0 galaxies, whose initial Q{r) is shown 
by a black solid line. Different orientation between the orbital 
angular momentum and the galaxy spins have been taken into 
account in this case (direct and retograde orbits), as well as dif- 
ferent orbital energies and angular momenta (orbit 01 and 05). 



Appendix A: Mixing processes and their 
dipendence on the orbital parameters 



In Sect |4.1| we have seen that the evolution of the metallicity 
profile for the fiducial merger between two gEO galaxies does 
not depend on the energy and angular momentum of the orbit. 
As discussed in Sect.|2l the gEO models represent quite idealized 




^ final 

Fig. A.3. Redistribution of stellar particles in the final remnants 
of gEO-gEO mergers. For every stellar particle in the remnant, 
we have evaluated its position in the remnant (through its dis- 
tance r final from the remnant center) and the position it had in 
the progenitor galaxy (through its distance r,„„ from the progen- 
itor center). The averaged quantities are then plotted, each curve 
corresponding to a different orbit, as explained in the legend of 
Fig lA.2| left panel. The corresponding standard deviations are 
also shown. Note that particles are redistributed in a very similar 
way in the merger remnants (similar averages and dispersions), 
independent of the orbit of the encounter. 



initial conditions, being spherical systems, without rotation. The 
final product of a merger of two not rotating spherical galaxies 
is typically a triaxial system, with some amount of rotation, due 
to the redistribution of the orbital angular momentum in internal 
one. Thus it can be interesting to understand if "hierarchical" 
initial conditions (the initial galaxies being the remnants of 
previous mergers) still produce final metallicity profiles which 
show similar properties, independently on the parameters of the 
interaction. A different alignement between the initial galaxy 
spins and the orbital angular momentum can lead, for example, 
to different tidal effects in the interacting pairs. Can this also 
lead to final metallicity profiles depending on the initial orbital 
energy and angular momentum? In other words, are the results 
found in Sect |4. ll confirmed when repeated hierarchical mergers 
are considered? 

To assess the generality of the results, we first simulated the 
merger of two iEO galaxies (not rotating spherical galaxies 
whose masses are equal to half of the gEO mass) and then two 
identical copies of the corresponding remnants were merged 
again, changing the initial orbital energy and angular momen- 
tum and its orientation with respect to the galaxy spins. 
The result of these experiments are shown in Fig lA.H where 
the metallicity profiles of 2iE0-2iE0 mergers are shown, for 
different orbits, corresponding to different initial values of the 
orbital angular momentum (orbit id=01 and 05) and to different 
alignments of this orbital angular momentum with respect to the 
galaxy spins (dirOO and retOO corresponding to the most extreme 
cases, when the orbital angular momentum is respectively paral- 
lel or antiparallel to the galaxy spins). In all the cases examined, 
it results that the final metallicity profiles are independent on 
the orbit of the encounter, at least inside the half mass radius of 
the remnant galaxy. It is interesting to note that this conclusion, 
i.e. the independence of the remnant metallicity profiles on the 
orbit of the encounter, does not depend on the morphological 
and kinematical properties of the initial systems either Indeed 
it is valid when the merger involves two spherical symmetric, 
not rotating systems (Sect |4.1| ), as well as when it involves two 
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triaxial, rotating galaxies (Fig lA.U . 

But where does this behaviour come from? In other words, 
how do mixing and phase-space density evolution proceed in 
collisionless major mergers and how do they depend on the or- 
bital parameters of the encounter? The collisionless Boltzmann 
equation tells us that the evolution of a collisionless system is 
characterized by the conservation of the six-dimensional phase- 
space density (also known as the "fine-grained" phase-space 
density). The computation of this phase-space density being 
quite difficult, one usually computes the "coarse-grained" phase- 
space density, which represents the average of the fine-grained 
one over some finite volume in the phase-space, a nd which obeys 
one o f the Mixing Theorems (see for example iTremaine et alj 
Il986l) or the spherically averaged quantity Q(r) - p(r)/cr'(r) - 
being p{r) the spherically averaged volume density and cr(r) the 
velocity dispersion- which still has the dimensions of a phase- 
space density, but being only a proxy of the coarse-grained distri- 
bution fun c tion, it does not satisfy any of the Mixing Theorems. 
IVass et all ( |2008|) have recently studied the evolution of the 
phase-space density distribution of dark matter halos during 
equal-mass mergers. They have taken into account a variety of 
shapes for the dark matter halos (from steep cusps to core-like 
profiles), of different orbital parameters, as well as recursive 
mergers of dark halos, in order to study the effect of a hierar- 
chical building up of halos on the evolution of the phase-space 
density profiles. They studied both the evolution of the coarse- 
grained distribution function F{r) and of its spherically averaged 
proxy Q{r), showing that they evolve in a similar way during the 
merging process, independently on the orbital properties of the 
encounter and on the number of repeated collisions. In all cases, 
the inner slopes and the overall shapes of the phase-space den- 
sity distribution of the merger remnant are close to that of the 
progenitor system. 

We too found similar results evaluating the Q{r) profile of 
our remnant galaxies, for different orbital initial conditions. 
As shown in Fig |A.2| mergers of spherically symmetric (upper 
panel) and of triaxial (lower panel) systems lead to an evolution 
of the spherically averaged Q{r) which is independent on the or- 
bit and which retains memory of the progenitor profile. 
This independency in the evolution of the phase-space distribu- 
tion on the orbital properties of the encounter leads also to very 
similar evolutions in the spatial distribution of particles after 
the collisions. In other words, stellar particles redistribute them- 
selves in the remnant galaxy in a very similar way, independently 
on the orbits (see Fig |A.3t , at least inside the remnant half-mass 
radius (rso - 7kpc for gEO-gEO mergers). This ultimately gives 
rise to the very similar metallicity profiles found in Fig. |2] and 
Fig.lAll 
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Table 4. Initial and final metallicity gradients for the fiducial 1:1 dry mergers. 
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Run Type 


Initial slope 
(Galaxy 1) 


Initial slope 
(Galaxy 2) 


Final slope 
(Remnant) 


Run Type 


Initial slope 
(Galaxy 1) 


Initial slope 
(Galaxy 2) 


Final slope 
(Remnant) 


gradOlgradOl 


-0.26 


-0.26 


-0.17 


grad05gradll 


0.00 


-0.30 


-0.11 


grad01grad02 


-0.26 


0.00 


-0.09 


grad05gradl2 


0.00 


-0.33 


-0.12 


grad01grad03 


-0.26 


-0.13 


-0.13 


grad05gradl3 


0.00 


-0.37 


-0.13 


grad01grad04 


-0.26 


-0.13 


-0.13 


grad05gradl4 


0.00 


-0.41 


-0.13 


grad01grad05 


-0.26 


0.00 


-0.08 


grad05gradl5 


0.00 


-0.44 


-0.14 


grad01grad06 


-0.26 


-0.22 


-0.16 


grad06grad06 


-0.22 


-0.22 


-0.15 


gradOlgradO? 


-0.26 


-0.18 


-0.15 


grad06grad07 


-0.22 


-0.18 


-0.14 


gradOlgradOS 


-0.26 


-0.15 


-0.14 


grad06grad08 


-0.22 


-0.15 


-0.13 


grad01grad09 


-0.26 


-0.07 


-0.11 


grad06grad09 


-0.22 


-0.07 


-0.10 


gradOlgradlO 


-0.26 


-0.04 


-0.09 


grad06gradl0 


-0.22 


-0.04 


-O.OS 


gradOlgradl 1 


-0.26 


-0.30 


-0.18 


grad06gradl 1 


-0.22 


-0.30 


-0.17 


grad01gradl2 


-0.26 


-0.33 


-0.19 


grad06gradl2 


-0.22 


-0.33 


-O.IS 


gradOlgradl 3 


-0.26 


-0.37 


-0.20 


grad06gradl3 


-0.22 


-0.37 


-O.IS 


gradOlgradl 4 


-0.26 


-0.41 


-0.20 


grad06gradl4 


-0.22 


-0.41 


-0.19 


gradOlgradl 5 


-0.26 


-0.44 


-0.21 


grad06gradl5 


-0.22 


-0.44 


-0.19 


grad02grad02 


0.00 


0.00 


0.00 


grad07grad07 


-0.19 


-0.18 


-0.13 


grad02grad03 


0.00 


-0.13 


-0.08 


grad07grad08 


-0.19 


-0.15 


-0.11 


grad02grad04 


0.00 


-0.13 


-0.05 


grad07grad09 


-0.19 


-0.07 


-0.09 


grad02grad05 


0.00 


0.00 


0.00 


grad07gradl0 


-0.19 


-0.04 


-0.07 


grad02grad06 


0.00 


-0.22 


-0.12 


grad07gradl 1 


-0.19 


-0.30 


-0.15 


grad02grad07 


0.00 


-0.18 


-0.11 


grad07gradl2 


-0.19 


-0.33 


-0.16 


grad02grad08 


0.00 


-0.15 


-0.09 


grad07gradl3 


-0.19 


-0.37 


-0.17 


grad02grad09 


0.00 


-0.07 


-0.06 


grad07gradl4 


-0.19 


-0.41 


-0.17 


grad02gradl0 


0.00 


-0.04 


-0.03 


grad07gradl5 


-0.19 


-0.44 


-O.IS 


grad02gradl 1 


0.00 


-0.30 


-0.14 


grad08grad08 


-0.15 


-0.15 


-0.10 


grad02gradl2 


0.00 


-0.33 


-0.15 


grad08grad09 


-0.15 


-0.07 


-0.08 


grad02gradl3 


0.00 


-0.37 


-0.16 


gradOSgradlO 


-0.15 


-0.04 


-0.06 


grad02gradl4 


0.00 


-0.41 


-0.17 


gradOSgradl 1 


-0.15 


-0.30 


-0.14 


grad02gradl5 


0.00 


-0.44 


-0.17 


grad08gradl2 


-0.15 


-0.33 


-0.14 


grad03grad03 


-0.13 


-0.13 


-0.09 


grad08gradl3 


-0.15 


-0.37 


-0.15 


grad03grad04 


-0.13 


-0.13 


-0.07 


grad08gradl4 


-0.15 


-0.41 


-0.16 


grad03grad05 


-0.13 


0.00 


-0.04 


gradOSgradl 5 


-0.15 


-0.44 


-0.16 


grad03grad06 


-0.13 


-0.22 


-0.11 


grad09grad09 


-0.07 


-0.07 


-0.05 


grad03grad07 


-0.13 


-0.18 


-0.11 


grad09gradl0 


-0.07 


-0.04 


-0.04 


grad03grad08 


-0.13 


-0.15 


-0.09 


grad02gradl 1 


-0.00 


-0.30 


-0.14 


grad03grad09 


-0.13 


-0.07 


-0.07 


grad09gradl2 


-0.07 


-0.33 


-0.11 


grad03gradlO 


-0.13 


-0.04 


-0.06 


grad09gradl3 


-0.07 


-0.37 


-0.11 


grad03gradl 1 


-0.13 


-0.30 


-0.13 


grad09gradl4 


-0.07 


-0.41 


-0.12 


grad03gradl2 


-0.13 


-0.33 


-0.14 


grad09gradl5 


-0.07 


-0.44 


-0.12 


grad03gradl3 


-0.13 


-0.37 


-0.14 


gradlOgradlO 


-0.04 


-0.04 


-0.03 


grad03gradl4 


-0.13 


-0.41 


-0.15 


gradlOgradll 


-0.04 


-0.30 


-O.OS 


grad03gradl5 


-0.13 


-0.44 


-0.15 


gradl0gradl2 


-0.04 


-0.33 


-0.09 


grad04grad04 


-0.13 


-0.13 


-0.09 


gradl0gradl3 


-0.04 


-0.37 


-0.09 


grad04grad05 


-0.13 


0.00 


-0.04 


gradl0gradl4 


-0.04 


-0.41 


-0.10 


grad04grad06 


-0.13 


-0.22 


-0.14 


gradl0gradl5 


-0.04 


-0.44 


-0.10 


grad04grad07 


-0.13 


-0.18 


-0.13 


gradllgradll 


-0.30 


-0.30 


-0.19 


grad04grad08 


-0.13 


-0.15 


-0.11 


gradl lgradl2 


-0.30 


-0.33 


-0.20 


grad04grad09 


-0.13 


-0.07 


-0.08 


gradl 1 gradl 3 


-0.30 


-0.37 


-0.21 


grad04gradl0 


-0.13 


-0.04 


-0.06 


gradl 1 gradl 4 


-0.30 


-0.41 


-0.22 


grad04gradl 1 


-0.13 


-0.30 


-0.16 


gradl 1 gradl 5 


-0.30 


-0.44 


-0.22 


grad04gradl2 


-0.13 


-0.33 


-0.17 


gradl2gradl2 


-0.33 


-0.33 


-0.21 


grad04gradl3 


-0.13 


-0.37 


-0.18 


gradl2gradl3 


-0.33 


-0.37 


-0.22 


grad04gradl4 


-0.13 


-0.41 


-0.19 


gradl2gradl4 


-0.33 


-0.41 


-0.23 


grad04gradl5 


-0.13 


-0.44 


-0.20 


gradl2gradl5 


-0.33 


-0.44 


-0.24 


grad05grad05 


0.00 


0.00 


0.00 


gradl 3 gradl 3 


-0.37 


-0.37 


-0.23 


grad05grad06 


0.00 


-0.22 


-0.10 


gradl 3 gradl 4 


-0.37 


-0.41 


-0.24 


grad05grad07 


0.00 


-0.18 


-0.09 


gradl 3 gradl 5 


-0.37 


-0.44 


-0.25 


grad05grad08 


0.00 


-0.15 


-0.07 


gradl4gradl4 


-0.41 


-0.41 


-0.25 


grad05grad09 


0.00 


-0.07 


-0.05 


gradl4gradl5 


-0.41 


-0.44 


-0.26 


gradOSgradlO 


0.00 


-0.04 


-0.03 


gradl5gradl5 


-0.44 


-0.44 


-0.27 



